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1. Introduction 

Freeze out (FO) is a term referring to the stage of expanding or exploding matter 
when its constituents (particles) lose contact, collisions cease, and local dynamical 
equilibrium is not maintained. In local equilibrium the evolution of the system can 
be described by hydrodynamics, while as time passes the system becomes more 
dilute, the number of non-interacting particles increases until the whole system 
breaks up and the particle momentum is frozen out. 

2. Freeze out and the Boltzmann Transport Equation 

The conservation of the total number of particles with collisions among them re- 
quires the relativistic Boltzmann Transport Equation (BTE), 

P^f = \f vffldihW& -\J 2V 34 ff 2 W$ , (1) 

commonly written in form of gain and loss terms. We have used the following short- 
hand notation for the invariant scalar product a bT^ c = ^t^- Furthermore, 
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the dynamics is governed by the invariant transition rates, Wjji, which stands for 
the elementary reaction, a + b — > c + d, satisfying the energy momentum conserva- 
tion, while fj = f(x,pj). 

To describe a gradual freeze out process, we split the / distribution function 
into two parts: / = f + f? . The result of collisions is the drain of particles from 
the interacting component, /*', which gradually builds up the free component, 
expressed by the freeze out probability, Vf. The FO probability populates the free 
component, while the rest, (1 — Vf), feeds the interacting component. We can 
separate the two components into two equations: 

jfW =+\j PfWg , (2) 

p^nf = -\J v2D i fiftp f wg+\ J uvjifiwg \ J 2 -D 3i rm$, 

so that the sum of these two equations returns eq. QJ. The free component does 
not contain a loss term because the free component cannot loose particles due to 
collisions. The latter two terms in the interacting component do not include Vf as 
they influence only the interacting term by redistributing particles in the momentum 
space and driving the interacting component towards rethermalization [[I]. 



3. Approximate kinetic freeze out equations 

We do not possess enough information to calculate the gain term containing the 
FO probability from eq. thus we rather approximate this term based on fun- 
damental physical principles with the so called escape rate, V eS c- The escape rate 
includes Vf, and separates the outgoing (gain) particles into a fraction that is still 
colliding and a fraction that is not. 

The probability not-to-collide with anything on the way out, should depend on 
the number of particles, which are in the way of a particle moving outwards in the 
direction p x /\p\ across the FO layer of thickness L. Following a particle moving 
outwards form the beginning, i.e. (x M = 0), to a position x**, with momentum p^, 
the actual distance up to the outer FO layer is (L — s)/ cos 0p = (L — x^da^)/ cos 6p, 
where s = x^da^, see Fig. ^ Assuming that the FO probability is inversely pro- 
portional to the remaining distance, the escape rate is: 

where the cut-off factor, Q(p fl d<r fl ), eliminates the FO of particles with negative 
momentum and A is the initial characteristic length of the system. If the particle 
momentum is not normal to the surface, the remaining spatial distance further 
increases with increasing 8^ and the probability to leave the system decreases as 
cos 6p, where 6p is the angle between the FO normal vector and p, see Fig. Q] This 
"naive" rate, cos Op, is not covariant, and did not take into account that the escape 
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rate of particles is also proportional to the particle velocity, \v\. The covariant 
generalization can be given in the following way 

cos^xN^^xM^f^V (4 ) 

where p^ = (p°,p) is the particle 4-momentum , u M is the flow velocity. Now, using 
the invariant form of the escape rate, the approximate equations governing the FO 
development for both time-like and space-like FO situations are: 

*"W - 4(rr|^:)(^f>^'+W.-/']i, 

where the equations describe a space-like, da^da^ = —1, or a time-like, da^da^ = 1, 
FO process depending on the FO normal vector. The equilibrium distribution func- 
tion is denoted by /* , while Xth is the rethermalization length or time. 




Fig. 1. Left: The picture of a gradual FO process within the finite FO layer, in 
coordinate space. The inside boundary of the FO layer, Si, indicates the points 
where the FO starts. This is the origin of the coordinate vector, x M . Within 
the finite thickness of the FO layer, L, the density of the interacting particles 
gradually decreases and disappears at the outside boundary S2 of the FO layer. 
Right: A simple hypersurface in the Rest Frame of the Gas (RFG: [t,x]) , where 
= (1,0,0,0), including time-like and space-like parts. On these two parts of 
the hypersurface in the Rest Frame of the Front (RFF: [t',x']), the normal to the 
hypersurface, da lll points into the direction of the t' (x')-axis respectively. 

We can study this new angular factor, W(p) — p^da^/p^UaQijp^da^), by taking 
different typical points of the FO hypersurface. At points A, B, C, the hypersurface 
is time-like, while at points D, E, F, the hypersurface is space-like. The resulting 
phase-space escape rates are shown in Fig. Elfor the cases B, C, D, E, F, in the 
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Rest Frame of the Gas (RFG) as well as in the Rest Frame of the Front (RFF). For 
point A on the hypersurface, W(p) = 1, in both reference frames. 

The effect of this relativistically invariant angular factor leads to a smoothly 
changing behavior of W(p), as the direction of the normal vector changes in RFG. 
In RFF, W(p) behaves discontinuously when we cross the light cone, from point C 
to point D, (see the contour line belonging to W(p) = 0). This is a consequence the 
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Fig. 2. The contour plots of the momentum dependent part of the escape rate, 
W(p), presented at different points of the FO hypersurface. The upper plots show 
the values of W(p) in RFG, while the plots below are in RFF. The momenta are in 
units of particle mass. 

The aim of freeze out calculations is to find the post FO momentum distribution 
and the relevant quantities from the properties of the matter on the pre FO side. 
The final post FO particle distributions are non-equilibrated and anisotropic distri- 
butions. These distributions in general cannot be Lorentz transformed to a frame 
where, the distribution is isotropic [EH S]- The only exception is when the normal 
to the FO hypersurface is parallel to the local flow velocity. The usual practice of 
assuming the Juttner or the cut-Juttner distribution as the post FO distribution 
is generally not valid. By introducing a finite thickness FO layer, we are strongly 
affecting the evolution of both the interacting and frozen out components [El E]- 
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